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ABSTRACT 

This  thesis  presents  an  analysis  of  the  throughput  characteristics  for  several 
classes  of  Aloha  packet  networks.  Specifically,  the  throughput  for  variable  packet 
length  Aloha  utilizing  multiple  power  levels  to  induce  receiver  capture  is  derived. 
The  results  are  extended  to  an  analysis  of  a  selective-repeat  Automatic  Repeat 
Request  (ARQ)  Aloha  network.  Analytical  results  are  presented  which  indicate 
a  significant  increase  in  throughput  for  a  variable  packet  network  implementing 
a  random  two  power  level  capture  scheme.  Further  research  into  the  area  of  the 
near /far  effect  on  Aloha  networks  is  included.  Improvements  in  throughput  for 
mobile  radio  Aloha  networks  which  are  subject  to  the  near/far  effect  are  presented. 
Tactical  Command,  Control,  and  Communications  (C3)  systems  of  the  future  will 
rely  heavily  on  Aloha  ground  mobile  data  networks.  The  incorporation  of  power 
capture  and  the  near/far  effect  into  future  tactical  networks  will  result  in  improved 
system  analysis,  design,  and  performance. 
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I.  INTRODUCTION 

United  States  Department  of  Defense  (DOD)  Tactical  Command,  Control,  and 
Communications  (C3)  systems  of  the  future  will  rely  heavily  on  ground  mobile 
data  networks.  These  networks  must  provide  high  reliability  and  connectivity  in  an 
environment  in  which  enemy  counter-C3  assets  will  attempt  to  disrupt  and  deny 
friendly  C3  networks.  A  great  deal  of  interest  has  been  focused  on  Aloha  network 
architectures  for  future  mobile  tactical  data  systems  due  to  their  simplicity  and 
robustness. 

Research  on  Aloha  has  been  conducted  since  the  early  1970s  [Ref.  l:p.  253].  Most 
of  this  work,  however,  has  utilized  idealized  channel  models.  While  this  has  pro- 
duced a  wealth  of  information  about  the  general  characteristics  of  Aloha  networks, 
it  remains  critical  that  real  world  conditions  be  taken  into  account.  Factors  which 
affect  Aloha  channels  in  the  real  tactical  environment  include  fading.  Near/Far  ef- 
fect, power  capture  (both  natural  and  induced),  and  jamming.  This  thesis  presents 
analysis  of  Near/Far  effect  and  power  capture  on  channels  which  employ  the  Aloha 
nei  work  discipline. 

Chapters  II  and  III  concentrate  on  the  effects  of  induced  power  capture  when 
employed  as  part  of  a  variable  packet  Aloha  network.  Chapter  IV  presents  an  anal- 
ysis of  the  ramifications  of  Near/Far  effect  on  three  different  Aloha  configurations. 
This  work  is  geared  specifically  toward  ground  mobile  data  network  applications. 


II.  POWER  CAPTURE  ALOHA 

Aloha  packet  networks  have  been  utilized  for  some  time  in  data  communications 
systems  which  are  composed  of  many  uncoordinated  users,  each  with  low  data  rate 
requirements.  Recent  emphasis  on  military  and  commercial  ground  mobile  data 
networks  will  increase  the  number  of  systems  which  employ  variations  of  the  Aloha 
protocol  as  a  channel  discipline.  Although  many  Aloha  systems  have  been  effec- 
tively employed  to  date,  several  variations  on  the  basic  Aloha  scheme  still  remain 
unstudied.  The  preponderance  of  work  on  Aloha  network  analysis  has  focused  on  an 
idealized  theoretical  radio  channel,  thereby  neglecting  the  effect  of  real  world  chan- 
nel conditions  on  system  performance.  It  is  particularly  important  to  incorporate 
channel  characteristics  into  Aloha  analysis  because,  in  contradiction  to  intuition, 
the  overall  throughput  of  an  Aloha  channel  can  actually  increase  under  certain  con- 
trived or  natural  channel  conditions.  A  specific  example  of  this  is  the  case  of  power 
capture  Aloha. 

Another  area  of  neglect  in  terms  of  Aloha  network  analysis  is  that  of  variable 
packet  Aloha.  Although  model  results  have  shown  that  fixed  packet  and  slotted 
Aloha  both  exhibit  better  throughput  characteristics  than  variable  packet,  the  real- 
ity of  bit  stuffing  logic  in  fact  reduces  the  actual  throughput  of  these  disciplines.  In 
contrast,  variable  packet  Aloha  can  be  shown  to  be  a  simple  and  efficient  network 
discipline. 

This  chapter  introduces  and  develops  the  throughput  characteristics  of  a  variable 
packet  Aloha  network  which  incorporates  the  power  capture  phenomenon.  For  this 
analysis,  only  created  capture  is  discussed,  although  the  results  can  easily  be  adapted 
to  accommodate  natural  capture  due  to  fading. 


Created  capture  effects  occur  when  groups  of  users  are  assigned  different  transmit- 
ter power  levels  in  an  effort  to  create  priority  classes,  or  when  all  users  dynamically 
select  a  random  transmit  power  in  order  to  increase  channel  throughput.  Previ- 
ous analysis  of  power  capture  in  slotted  and  fixed  packet  networks  has  shown  a 
dramatic  increase  in  channel  throughput  over  the  idealized  channel  models  [Ref.  2]. 
This  chapter  investigates  the  effect  of  created  power  capture  with  two  random  power 
levels  on  an  Aloha  network  using  variable  length  packets. 

A.      SYSTEM  DYNAMICS 

In  this  analysis,  it  is  assumed  that  there  is  an  infinite  user  population  and  that 
the  channel  traffic  rate  is  Poisson  with  parameter  g  packets  per  second.  Assume 
newly  created  packets  are  of  variable  length  x  (bits)  and  are  independently  and 
identically  distributed  according  to  the  probability  density  function  a(x)  with  mean 
x.  Within  this  model,  all  channel  packets  are  indistinguishable.  That  is  to  say  that 
newly  arrived  packets  and  retransmitted  packets  have  their  length  x  redrawn  afresh 
from  a(x)  for  every  packet.  This  is  an  analytical  device  known  as  the  independence 
assumption. 

The  probability  that  n  packets  interfere  (overlap)  with  an  arbitrary  tagged  packet 
during  its  transmission  is  given  by 

Pr{n}=  Wx  +  uWt-9i*+u)  (2.1) 

n\ 

where  x  is  the  tagged  packet  length,  u  represents  the  random  length  of  a  preceding 
packet,  and  n  is  the  number  of  interfering  packets  [Ref.  3].  For  the  case  of  n  =  0, 
this  can  be  interpreted  as  the  probability  that  the  tagged  packet  encounters  no 
overlap  from  preceding  packets  and  the  next  interarrival  is  larger  than  x.  Given  the 


case  where  no  interferers  can  be  tolerated  by  a  tagged  packet,  equation  (2.1)  can  be 
viewed  as  the  probability  of  successful  tagged  packet  transmission. 

This  analysis,  however,  considers  the  case  when  a  receiver  can  be  captured  by  a 
tagged  packet  in  the  presence  of  interfering  packets  if  the  tagged  packet  signal  power 
exceeds  the  joint  signal  power  of  all  interfering  packets  by  some  specific  capture 
threshold,  70.  If  at  any  time  during  the  tagged  packet  interval,  the  ratio  of  tagged 
packet  power  to  joint  interfering  power  falls  below  70,  the  tagged  packet  is  considered 
to  be  destroyed. 

Unlike  slotted  Aloha  where  the  number  of  interferers  is  constant  over  the  entire 
tagged  packet  transmission  interval  (Figure  2.1(a)),  or  fixed  packet  Aloha  where  the 
number  of  interferers  is  a  stochastic  process  made  up  of  early  and  late  interfering 
packets  (Figure  2.1(b)),  the  number  of  interferers  present  in  a  variable  packet  Aloha 
channel  represents  a  stochastic  process  which  is  driven  by  four  classes  of  interferers. 
Figure  2.1(c)  depicts  the  four  possible  interferers  types.  The  terms  early  and  late 
refer  to  the  begin  and  end  time  relationship  between  an  interfere!-  and  the  tagged 
packet.  For  instance,  the  interferer  identified  as  Early-Early  in  Figure  2.1(c)  cor- 
responds to  any  interferer  which  begins  before  the  tagged  packet  (early)  and  ends 
before  the  tagged  packet  (early). 

The  dynamics  of  the  interfering  stochastic  process  is  as  follows.  We  note  that  the 
arrival  and  departure  of  n  interfering  packets  partitions  the  tagged  packet  interval 
into  several  non-  overlapping  intervals  of  random  length.  The  exact  number  of  these 
intervals,  given  a  known  number  of  interferers  and  their  corresponding  interferer 
type,  is  found  by  the  following  relationship. 

Total   Intervals   =        (number   of   Late-Late  packets) 

+    (number   of   Early-Early  packets) 
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Figure  2.1:   Aloha  Interfering  Packets 


+   2  x    (number  of   Late-Early  packets)    +    1 


Relating  this  equation  to  Figure  2.1(c),  it  is  evident,  for  instance,  that  a  single  Early- 
Early  or  Late-Late  packet  would  partition  the  tagged  packet  into  1  +  1  intervals.  One 
interval  occurs  during  the  absence  of  the  interferer;  the  other  during  the  presence 
of  an  interferer.  In  a  different  case,  a  single  Late-  Early  would  partition  the  tagged 
packet  into  2  -f  1  intervals.  Two  of  the  partitions  would  have  no  interferers  present ; 
the  other  partition  would  occur  during  the  period  of  the  interfering  packet.  The 
Early-Late  packet  covers  the  entire  duration  of  the  tagged  packet  and  therefore 
does  not  provide  any  partitioning. 

Each  interval  has  associated  with  it  a  total  number  of  interferers  /(/,•),  where  /,  is 
the  start  time  of  the  i'th  partition  during  the  tagged  packet  interval.  Each  unique 
sequence  of  total  interferers  per  partition  (/(/,)  for  all  ?)  shall  be  called  a  realization. 
Figure  2.2  depicts  a  specific  realization  of  interfering  packets  for  a  case  where  n  =  3. 


TAGGED 
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Figure  2.2:   An  Example  of  an  Interference  Realization 

For  this  specific  realization,  the  tagged  packet  total  intervals  is  5.  The  total  number 
(I(t,))  of  interferers  present  during  each  interval  is  indicated  at  the  bottom  of  the 
figure. 


B.      MAXIMUM  NUMBER  OF  INTERFERERS 

A  computer  program  (Appendix  A)  was  written  to  generate  a  list  of  every  possible 
interfering  realization  for  a  given  n.  From  this  list  of  all  possible  realizations,  the 
number  of  realizations  (Cj)  in  which  the  maximum  number  of  interferers  present 
at  any  time  during  the  tagged  packet  interval  does  not  exceed  some  number  j  was 
tallied.  Table  2.1  gives  the  results  for  n  =  1,2,3,4,5. 

An  example  of  how  one  would  generate  a  Cj(n)  value  is  as  follows.  Consider 
the  simple  case  of  n  =  1.  In  this  case,  only  four  realizations  are  possible,  namely 
the  interfering  packet  can  be  an  Early-Early,  Early-Late,  Late- Early,  or  Late-  Late. 
Continuing  with  this  simple  example,  it  can  be  seen  that  given  there  is  1   known 


n 

Ci{n) 

C2(n) 

C3(n) 

C4(n) 

CB(n) 

/t 

f(n) 

1 

4 

0 

0 

0 

0 

4 

0.25 

2 
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10 

0 

0 
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14 

0.0714 

3 

4 

24 

20 

0 

0 

48 

0.0208 

4 

4 

50 

75 

35 

0 

1G4 

0.0061 

5 

4 

100 

225 

176 

5G 

561 

0.0018 

TABLE  2.1:   Maximum  Number  of  Interfering  Packets  Given  n 


interfere!-,  there  is  no  realization  where  the  maximum  number  of  simultaneous  in- 
terferers  during  the  entire  duration  of  the  realization  is  less  than  or  equal  to  j  =  0. 
In  fact,  in  this  example,  the  maximum  number  of  simultaneous  interferers  for  each 
of  the  four  possible  realizations  is  1.  Therefore,  the  value  of  C\(n)  for  n  —  1  is  4 
as  seen  in  Table  2.1.  For  the  case  of  n  =  2,  there  are  14  combinations  of  4  packet 
types  taken  2  at  a  time.  These  14  pairs  of  interfering  packets  generate  14  unique 
realizations.  Of  these  14  realizations,  4  have  a  maximum  simultaneous  number  of 
interferers  of  1,  and  10  have  a  maximum  of  2  interferers.  The  program  in  Appendix 
A  is  used  to  automate  the  process  of  generating  combinations  of  packets,  construct- 
ing the  resulting  realizations,  and  determining  the  maximum  number  of  interferers 
per  realization.  Table  2.1  is  the  result  of  the  program. 

Each  Cj(n)  gives  the  number  of  realizations  in  which  the  maximum  number  of 
interferers  at  any  time  equals  j,  given  the  total  number  of  interfering  packets  is 


known  to  be  n.  The  function  It  is  simply  the  total  number  of  all  possible  realizations 
given  n  packets  are  known  to  interfere.  The  function  f(n)  is  the  inverse  of  Ij  and 
is  simply  used  as  a  shorthand  notation  for  I /It- 

C.     VARIABLE  PACKET  ALOHA  WITH  RANDOM  POWER  LEV- 
ELS 

Data  packets  arrive  at  the  receiver  after  being  spatially  attenuated.  Each  packet 
is  transmitted  with  one  of  two  normalized  power  levels  given  by  the  set  Q,  =  {1,  M}. 
Here,  M  is  some  integer  multiple  of  the  threshold  70  times  the  lower  power  level 
P  —  1  (normalized),  that  is  M  =  N~f0  where  TV  is  an  integer  greater  than  1.  Each 
user  selects  a  transmit  power  level  for  each  individual  packet  from  the  set  £1  with 
equal  probability  ^  in  order  to  avoid  any  class  prioritization  among  user  groups. 
The  higher  power  level  M  is  chosen  according  to  the  relation 

(Ar+l)7o>  M  >  W70  (2.2) 

where  A'  >  1  is  an  integer  and  70  is  the  power  capture  threshold  of  the  receiver. 

The  tagged  packet,  arriving  at  the  receiver  with  power  Pt  E  {1,M}  may  capture 
the  receiver  given  a  realization  of  n  interfering  packets  if  and  only  if 

i(u) 

Pt  >  lomax[J2  Pj]  (2.3) 

3=1 

for  all  intervals  (fs)  in  the  realization  where  P3  is  the  power  level  of  packet  j . 

For  a  tagged  packet  to  capture  the  receiver,  two  events  must  occur.  First,  the 
tagged  packet  must  be  of  power  M  while  none  of  the  interferers  are  of  power  M. 
Call  this  event  A,  then 


Pr{A\n]  =  — (2.4' 


8 


The  second  condition  for  capture  is  that  the  tagged  packet  must  satisfy  (2.3)  where 
Pt  —  M  and  all  Pj  =  1.  Given  a  specific  realization,  this  equates  to 

Ar  >  rnax{I{ti))  (2.5) 

The  value  N  is  therefore  the  maximum  number  of  interferers  that  can  be  tolerated 
at  any  instant  of  time  during  the  tagged  packets  transmission  interval.  If  event  B 
denotes  the  tagged  packet  power  satisfying  (2.5),  then 

Pr{5|n}  =  /(n)f;CJ-(n)  (2.6) 

where  C3{n)  is  the  number  of  realizations,  given  n  interfering  packets,  in  which  the 
maximum  number  of  simultaneous  interferers  equals  j.  Table  2.1  provides  a  list  of 
Cj(n)'s  and  the  inverse  /(??)  of  the  total  number  of  possible  realizations  given  n 
interfering  packets. 

By  using  equations  (2.4  )  and  (2.6)  the  conditional  capture  probability  Pr{capture\n) 
is  obtained. 


Pr{capture\n,n  =  0}     =     1  (2.7) 

Pr{capture\n,n>  1}     =     Pr{A\n}Pr{B\n}  (2.8) 

=     /(n)2-("+1>]TC»  (2.9) 

3=0 

The  overall  probability  of  receiver  capture  by  a  tagged  packet  is  found  by  mul- 
tiplying the  conditional  probability  in  equation  (2.9)  by  the  probability  that  there 
are  n  packets  in  the  channel,  and  summing  over  all  n.  Equation  (2.1)  provides  the 
probability  of  n  packets,  therefore,  the  probability  of  receiver  capture  by  a  tagged 
packet  which  shall  be  called  the  probability  of  success  is 


Pr{success}  =  c-'<*+«>[1  +  £  /(n)2-^+1)Ml+^  £  C.(n)]  (2.i0) 


n! 

n=l  j=0 


The  conditional  channel  departure  rate  in  the  interval  x  +  u  is  given  by  Reference  3 
to  be 

<f)(x,u)  =  g(x -\- u)Pr {success}  (2-11) 

The  channel  throughput  5,  in  packets  per  mean  packet  time,  is  related  to  the  channel 
departure  rate  by 

/•oo     /-oo  2 

S=  I      /      -(j>{x,u)a(x)a(u)dxdu  (2.12) 

Jo     Jo     2 

Define  G  to  be  the  attempted  packet  rate  in  packets  per  mean  packet  time  (i.e. 
G  —  gx).  Appendix  B  contains  the  detailed  solution  of  the  above  equation.  The 
resulting  channel  throughput  is  found  to  be  the  following: 


5  =  TgTTp  + 1  |^/(»)Cw,(e  +  ir"tI|Ec,W  (2.13) 

\  '        /  n—\  j=0 

Figure  2.3  shows  the  throughput  of  variable  packet  Aloha  with  random  two  power 
signal  for  values  of  N  =  0.1,2,3.  Packet  lengths  are  exponentially  distributed  and,  for 
computational  purposes,  n  is  taken  to  be  between  1  and  5.  The  curve  corresponding 
to  N  =  0  is  for  conventional  variable  packet  Aloha. 
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Figure  2.3:  Throughput  of  Variable  Packet  Aloha  With  Power  Capture 
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III.  SELECTIVE-REPEAT  ALOHA  WITH 

CAPTURE 

To  continue  the  study  of  variable  packet  Aloha,  the  results  of  the  previous  chapter 
are  now  applied  to  a  more  advanced  Aloha  network  which  incorporates  selective- 
repeat  as  a  data  link  control. 

The  specific  architecture  discussed  here  is  a  network  where  remote  terminals  com- 
municate with  a  hub  station  via  an  Aloha  channel  in  the  query-response  mode.  The 
feedback  hub-to-remote  link  is  an  asychronous  time-division-multiplexed  (ATDM) 
channel.  Each  remote  terminal  receives  the  same  data  stream  and  selects  applicable 
packets  by  searching  for  its  own  address. 

Assume  that  packet  length  of  newly  arrived  and  retransmitted  packets  in  the 
Aloha  channel  are  independent  and  identically  distributed  with  probability  density 
function  a(x)  and  with  mean  x.  In  this  model,  packets  are  indistinguishable.  Both 
newly  arrived  packets  and  retransmissions  can  be  pictured  as  having  their  length 
redrawn  afresh  from  a(x)  prior  to  each  transmission. 

Packets  in  the  ATDM  channel  are  independent  and  identically  distributed  with 
density  function  b{y)  and  mean  y. 

Variable  packets  on  the  Aloha  channel  are  broken  up  into  z  numbered  and  con- 
stant length  minipaekets  prior  to  transmission.  When  packet  damage  occurs,  only 
the  specific  minipaekets  which  are  unreadable  are  retransmitted.  Let  D(z),z  =1,2,.. 
be  the  distribution  of  the  minipaekets.  On  the  average,  each  large  packet  is  seg- 
mented into  z  minipaekets. 
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Let  the  channel  attempted  packet  rate  be  g  packets  per  second  and  zg  the  channel 
attempted  rate  in  minipackets  per  second.  With  m  being  the  fixed  length  of  a 
minipacket,  the  attempted  channel  rate  in  minipackets  per  minipacket  time  is  rnzg. 

The  conditional  probability  Ps  that  a  minipacket  encountered  no  collisions  is  the 
probability  of  no  overlap  from  preceding  packets  and  that  the  next  interarrival  is 
larger  than  the  minipacket  size  m.  Thus 

Pa(z)  =  e-9(mz+m)  (3.1) 

In  a  capture  environment,  however,  Ps  must  take  into  account  the  times  when  the 
tagged  packet  captures  the  receiver  in  spite  of  the  presence  of  n  interferers.  There- 
fore, using  equation  (2.1)  from  the  previous  chapter,  the  probability  of  minipacket 
success  is 

N 


Ps{z)  =  c-*(-+»)[]  +  t  finV-wMrnz  +  m)}"  £ 

n=l  j=0 

The  conditional  probability  of  successful  minipacket  transmission  is 


q(z,y)  =  Ps(z)[l-Pl(m)]ll-p2(y)}  (3.3) 

where  Pi(m)  is  the  minipacket  error  probability  due  to  errors  in  the  Aloha  channel 
and  P2{y)  is  the  error  rate  in  the  ATDM  channel.  The  conditional  minipacket 
departure  rate  in  the  interval  of  z  +  1  minipackets  is 

(t>m{z,y)  =  mg(z  +  l)q{z,y)  (3.4) 

The  average  minipacket  throughput  is  given  by  the  following  equation: 

z       f°° 

Sm  =  z— —  /     4>m{z,y)dy  (3.5) 

z  -\-  1  Jo 
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Assume  that  each  packet  consists  of  a  geometrically  distributed  number  of  fixed 
length  minipackets,  then  the  distribution  D(z),  which  is  the  probability  that  a  packet 
consists  of  z  minipackets,  is  defined  as  follows. 

/?(*)  =  (1  -  M)-V  (3.6) 

where  n  is  the  probability  that  a  rninipacket  is  generated.  Because  ~z  =      *   ,,  D(z) 
can  be  written  as 

D{z)  =  \{l-\)n-1  (3.7) 

Appendix  C  provides  the  detailed  solution  for  the  throughput  equation.  The 
resulting  rather  lengthy  set  of  equations  is  as  follows: 


Sm  =  A  +  ^~a(l  -  Pl(m))  £)  h(n)w-7[F(n)R(n)  -  8}  (3.8) 

~  +  *■  n=l 


4       oGv2{a-v){\  -pi{m)) 

z{z  +  l)w2(l-v)2  [       ' 


F(n)  =  (1  -  u)-<*+1)(i  +  1)!  (3.10) 


R(n)  =  J2vrT,(-l)kE(k,r)  (3.11) 

r  =  l  fc=0 


E{k,r)=  ,„        7    '  (3.12) 


ib!(t-Jfc  +  l)! 

G  =  ml5  (3.13) 
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fc(n)  =  /(„)2-("+1)-1^Ty  £  <?»  (3.14) 


N 

n\z<~ 

j=0 


u=  (l_I)e-(f)  (3.15) 

2 


u>  =  (l-z)  (3.16) 

z 


k  =  z  +  l  (3.17) 

<  =  n  +  l  (3.18) 

Figure  3.1  shows  the  throughput  for  the  variable  packet  Aloha  with  selective- 
repeat  data  link  control.  The  parameters  are:  p\(m)  =  0,  p^i™)  —  0,  x  and  y  are 
both  exponentially  distributed,  x  —  y ,  and  3  =  8  minipackets. 
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IV.  NEAR/FAR  EFFECT  ON  ALOHA 

MODELS 

This  final  analysis  chapter  deals  with  a  network  phenomenon  which  is  specific  to 
mobile  packet  radio  networks  that  do  not  incorporate  adaptive  control  of  terminal 
transmit  power.  The  near/far  effect  is  a  condition  which  occurs  when  packets  arrive 
at  a  receiver  with  different  mean  power  due  to  differences  in  transmission  distances. 
Under  these  conditions,  a  tagged  packet  may  be  able  to  capture  the  receiver,  in  spite 
of  interfering  packets,  if  the  tagged  packet  power  is  sufficiently  greater  than  the  joint 
power  of  the  interfering  packets.  Because  each  packet  is  spatially  attenuated  due 
to  the  near/far  effect,  all  packets  have  a  finite  probability  of  capturing  the  receiver. 
The  throughput  characteristics  of  various  Aloha  network  configurations  which  are 
subject  to  the  near/far  effect  are  derived  and  presented. 

A.      CELL  MODEL 

The  mean  power  w  of  a  packet  at  a  distance  7-  from  the  transmitter  is  of  the 
general  form 

w  =  cr~a  (4.1) 

In  the  event  of  ground-wave  propagation  without  shadowing. 

c  =  PTlGTlGRH2TlH2R  (4.2) 

where  Pjn  Gji-  and  Hjt  ai'e  the  transmit  power,  antenna  gain,  and  antenna  height 
(above  ground),  respectively,  of  the  transmitting  mobile  terminal.  The  values  Gr 
and  Hr  are  the  gain  and  height  above  ground  of  the  base  station  antenna.    The 


exponent  a  gives  the  attenuation  law  for  the  channel  considered  (2  <  a  <  5).  In 
the  event  of  UHF  propagation,  a  typical  value  is  a  =  4  [Ref.  3:p.  264]. 

Assuming  identical  mobile  terminals  and  omnidirectional  antennas  c  is  normal- 
ized to  unity  for  all  users  without  loss  of  generality  since  receiver  capture  power  is 
determined  by  the  ratio  of  signal  powers  [Ref.  3:p.  264].  Additionally,  the  value  of 
a  in  this  analysis  is  taken  as  4. 

In  this  model,  signal  power  depends  on  the  spatial  distribution  of  users  in  the 
network.  Users  create  a  traffic  density  g(r)  within  a  circular  cell  around  a  base 
station  receiving  node.  The  density  g(r)  can  be  thought  of  as  the  normalized  packet 
traffic  per  unit  area  at  a  distance  r  from  the  base  receiver.  The  packet  generation 
rate  in  a  given  area  depends  on  the  distance  to  the  base  station  r  and  is  independent 
of  direction.  Given  this  cell  model,  the  total  traffic  rate  offered  in  the  network  can 
be  described  as 

Jroo 
1     xg{x)dx  (4.3) 

o 

Figure  4.1  depicts  the  general  cell  model  used  in  this  analysis. 

The  spatial  distribution  of  packet  generation  is  found  to  be  the  probability  a 
packet  is  generated  within  a  distance  R.  That  is  to  say 

FR(r)  =  Pr{R  <  r)  =  ^  f  xg(x)dx  (4.4) 

G  Jo 

Differentiating,  the  probability  density  function  for  packet  generation  is  found  to 
be 


2tt 
fR(r)  =  -^rg(r)  (4.5) 
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Figure  4.1:   Cell  Model  for  Near/Far  Effect  Analysis 

As  a  means  of  simplification,  users  are  assumed  to  be  uniformly  distributed  in  a 
circle  of  radius  rmax  with  I  he  base  receiving  station  at  the  center  of  the  cell  as  shown 
in  Figure  4.2.  As  seen  in  Figure  4.2,  g(r)  is  normalized  to  G  l^r^nax  for  0  <  7-  <  rmax. 

From  this,  the  equations  for  the  total  network  traffic  and  packet  generation  reduce 
to 


G  =  TTT 


2 
»»i  ax 


(4.6) 


h(r)  = 


27rr 
~G 


(4.7) 


for  0  <  r  <  rmax 

A  standard  change  of  variable  operation  is  performed  on  equations  (4.1)  and  (4.6) 
to  define  a  probability  density  function  for  the  packet  power  iv. 
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Figure  4.2:  Distribution  of  Users  Within  the  Cell 


2tt 


fw{w)  =  —-w  (°  +  1) 

ah 


(4.8) 


for  rm°    <  w  <  oo. 


At  this  point,  the  network  cell  model  has  been  described  with  the  distribution 
of  packet  generation  given  by  equation  (4.6)  and  the  packet  power  density  function 
given  by  equation  (4.7).  The  near/far  effect  is  now  incorporated  into  the  model. 

B.      NEAR/FAR  EFFECT 

In  a  phenomenon  which  is  similar  in  many  respects  to  the  power  capture  case 
described  earlier,  the  near/far  capture  occurs  when  the  ratio  of  a  tagged  packets 
power  to  the  joint  interfering  power  exceeds  some  receiver  capture  threshold.  Again, 
the  number  of  interferers  during  the  tagged  packets  transmission  is  a  stochastic 
process.  The  capture  ratio  is  described  as 
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Z=y  (4.9) 

where  X  is  the  packet  power  of  the  tagged  packet  and  Y  is  the  joint  interfering 
power  of  some  number  n  interfering  packets.  As  in  previous  discussions,  the  receiver 
is  captured  if  the  ratio  described  by  Z  exceeds  the  threshold  given  by  the  parameter 

7o- 

The  distribution  of  X  is  simply  the  single  packet  power  density  function  given 
in  (4.7).  The  joint  interfering  power  density  Y  is  dependent  on  the  number  of 
interferers  and  is  given  as 

fY(y\n)  =  [fw{w)Y®  (4.10) 

where  n .;__,  denotes  the  ??-fold  convolution.  For  the  cases  of  it  =  0  and  n  —  1,  no 
convolutions  are  required  to  determine  fy{y).  For  higher  numbers  of  interferers,  a 
computer  based  convolution  package  is  used  to  generate  specific  values  for  fy(y\n). 
Appendix  D  provides  details  on  the  analytical  result  for  the  case  of  n  =  1  and  the 
numerical  results  for  7/  =  2.3. 

Utilizing  the  density  functions  for  A'  and  Y,  the  distribution  of  Z  can  be  found 
using  the  following  formula. 

fz{z\n)=  f°°  yfx(yz)fy(y\n)dy  (4.11) 

Jo 

The  capture  probability,  that  is  the  probability  that  Z  exceeds  the  receiver  thresh- 
old 70  given  n  interferers,  is  found  by  integrating  fz{z\n)  from  70  to  00. 

Table  4.1  presents  the  probability  of  capture  given  n  interferers  for  the  range 
0  <  n  <  4.  Although  the  case  of  ??  =  1  is  in  a  closed  form  with  parameter  70,  all 
subsequent  capture  probabilities  assume  a  value  of  70  =  3. 
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n 

Pr{capiurc\n} 

0 
J 
2 

1 

•2/0 

0.057 

3 

0.022 

4 

0.0J0 

TABLE  4.1:  Capture  Probabilities 

The  limited  list  of  capture  probabilities  listed  in  Table  4.1  is  sufficient  to  provide 
excellent  resolution  of  the  throughput  characteristics  for  the  Aloha  models  discussed 
in  the  following  sections.  The  probability  of  capture  for  cases  where  there  are  greater 
than  4  interfering  packets  is  so  small  as  to  be  considered  insignificant. 

C.      ALOHA  NETWORK  MODELS 

Three  common  Aloha  models  are  analyzed  with  respect  to  the  improvement  in 
overall  network  throughput  gained  from  the  near/far  effect.  In  each  case,  the  general 
definition  of  throughput  is 


S  =  GPr  {success} 
which  in  a  capture  environment  can  be  viewed  as 


(4.12) 
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S  =  GY,  Pr{capture\n}Pr{n]  (4.13) 

The  following  sections  use  this  approach  to  determine  the  network  throughput  for 
slotted  and  fixed  packet  Aloha  networks.  Variable  packet  Aloha  is  analyzed  in  the 
same  general  manner  as  was  done  in  Chapter  I  for  power  capture. 

1.  Slotted  Aloha 

For  slotted  Aloha,  the  probability  that  n  packets  are  generated  during  a 
given  packet  time  is  given  by  the  Poisson  Distribution  [Ref.  l:p.  255]: 

Pr{n)  =  ^-  (4.14) 

Therefore,  the  general  form  for  the  throughput  in  a  slotted  Aloha  network  with  the 
near/far  effect  is 

x                                Gne~G 
S  =  GY,  Pr{capture\n] —  (4.15) 

71=0  n- 

Using  the  Pr{capturt\n}  values  in  Table  4.1,  the  throughput  S  is  calculated  and  is 
shown  in  Figure  4.3. 

2.  Fixed  Packet  Aloha 

The  number  of  packets  generated  during  a  given  packet  time  in  a  fixed 
packet  network  is  given  by  [Kef.  4:p.   11]: 


Pr{n]  =  {^Le~2G  (4.16) 


resulting  in  a  throughput  of 


00  (2G)n 

S  =  G  J2  pr{capture\n}- — y-e~2G  (4.17) 

n  =  0  n- 
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Figure  4.3:  Slotted  Aloha  With  Near/Far  Effect 
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Figure  4.4:   Fixed  Packet  Aloha  With  Near/Far  Effect 

Again,  using  the  numerical  results  for  the  probability  of  capture  given  in  'Table 
4.1,  the  throughput  is  calculated  and  is  depicted  in  Figure  4.4. 
3.      Variable  Packet  Aloha 

This  analysis  is  very  similar  to  the  work  presented  in  Chapter  II  and 
therefore  the  variable  names  used  here  are  the  same.  First,  the  probability  of  success 
for  a  tagged  packet  is  defined  as 


Ps{x,u)  =  T  Pr{capiure\n}e-g{T+u)  ^  +  ^" 

n! 


n=0 


(4.18) 
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For  this  analysis,  perfect  channel  conditions  are  assumed  in  order  to  simplify  the 
calculations.  The  conditional  probability  of  a  successful  transmission  in  this  case  is 
seen  to  be 

q{x,u)  =  Ps(x,u)  (4.19) 

The  conditional  channel  departure  rate  in  the  interval  x  +  u  is 

<j){x,u)  =  g(x  -f  u)q(x  +  u)  (4.20) 

From  this,  the  channel  throughput  can  be  calculated  from  the  following  equation: 

fix      roo 

$  —  T  /       /      4>{x,u)a(x)a(u)dxdu  (4-21) 

2  Jo     Jo 

Take  the  case  where  a(x)  and  a(u)  are  exponentially  distributed,  the  above  equa- 
tion simplifies  to 

S  =  — ^—  +  £  '^llpricaptureWG^iG  +  1)'("+3)  (4.22) 

(6  + I)3      „=i        2 

Using  the  numerical  results  for  Pr{caTpture\n]  given  in  Table  2.1,  the  throughput 
for  variable  packet  Aloha  with  Near/Far  effect  is  shown  in  Figure  4.5. 
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Figure  4.5:  Throughput  for  Variable  Packet  Aloha  With  Near/Far  Effect 
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V.  CONCLUSIONS 

This  thesis  has  demonstrated  the  effects  of  induced  two  level  power  capture  on 
variable  packet  Aloha  networks  and  the  effects  of  the  near/far  phenomenon  on  var- 
ious Aloha  configurations.  In  all  cases,  significant  improvements  in  overall  network 
throughput  were  achieved. 

The  development  of  highly  reliable  Tactical  data  systems  of  the  future  depend 
on  Aloha  as  a  multiaccess  scheme  for  ground  mobile  terminals.  It  is  critical  that 
these  networks  be  allowed  to  achieve  their  maximum  potential.  Given  the  increased 
reliance  on  tactical  data  communications,  Aloha  networks  which  provide  the  maxi- 
mum user  throughput  possible,  even  in  a  jamming  environment,  must  be  developed. 

The  increases  in  overall  network  throughput  presented  in  this  thesis  indicated 
the  improvements  in  performance  which  can  be  realized  when  power  capture  and 
the  near/far  effect  are  incorporated.  As  seen  in  Chapter  2,  induced  power  capture 
increased  the  network  throughput  of  a  variable  packet  Aloha  network  from  0.148 
packets  per  mean  packet  time  in  the  idealized  case  to  0.184  in  a  capture  environment. 
Utilizing  the  selective-repeat  data  link  control  with  power  capture  variable  packet 
Aloha  increased  throughput  from  a  maximum  of  0.256  to  about  0.354.  In  the  final 
analysis  chapter,  improvements  in  throughput  for  slotted,  fixed  packet  and  variable 
packet  Aloha  which  incorporated  the  near /far  effect  were  presented. 

Continued  work  in  the  area  of  Aloha  network  analysis  will  be  needed  in  order 
to  develop  the  high  capacity,  high  reliability  tactical  data  networks  required  in  the 
future.  Incorporating  power  capture  and  the  near/far  effect  into  future  designs  will 
result  in  better  systems  analysis  and  design  for  these  tactical  mobile  networks. 
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APPENDIX  A 


{  } 

{  PROGRAM  CJ  } 

{  } 

{   Joe  McCartin,  January  1989  } 

{  } 

{  Turbo  Pascal  -  IBM  PC  } 

{  } 

{  Description:  } 

{  } 

{  This  program  is  designed  to  generate  all  } 

{  realizations  of  interferers  given  n  packets  are  } 

{  known  to  interfere.   Using  the  realizations,  a  } 

{  table  of  Cj (n)  values  is  compiled.   This  program  } 

{  is  geared  specifically  toward  the  purposes  of  this  } 

{  thesis  and  is  not  a  general  purpose  program.  } 

{  } 

-T  *  *  *  *  *  *  *  4  4  4  4  4  4  4  4  4    4  4   +  4  4  4   +  +  4  4   4  4  4   4    +  t  *  *  •  4  4   4   4  4  4   *  ♦  +  *  +  4  4  4   4   4  4   4   +  \ 

program  cj  (filevar,  output); 
const 

max_mtervals   =    11; 
matnx.length   =    1024; 
cj _matrix_dim   =    10; 

type 

realization.vector   =   array    [1 .  .max_mtervals]    of 

integer; 
realization_matrix   =   array    [1 .  .matrix_length]    of 

realization. vector; 
cj_vector  =   array    [1 . . cj_matrix_dim]    of    integer; 
cj .array   =   array    [l . . cj_matrix_dim]    of   cj .vector; 
index   =   array    [1..4]    of    integer; 
comb.matnx   =   array    [1..220]    of    index; 

var 

realization.x , 

realizat ion_y : realization.matnx ; 

vector : comb_matr ix ; 
vector.length : integer ; 

active_matrix , 
active_matrix_length: integer; 

c j .matrix: cj .array ; 
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init.interf erers , 
max_interf erers , 
min_interf erers , 
num_intervals , 
final. interf erers : integer ; 

n, 
i, 

J. 
t, 
h, 
runs : integer; 

file_name: string [10] ; 
f ilevar : text ; 
outf ilevar : text ; 

ans :char ; 

archive: boolean; 

{  } 

{  INITIALIZE.CJ  MATRIX                              } 

{  } 

{  Tabulated  values  of  Cj ' s  are  kept  in  a  square  matrix  } 

{  defined  by  the  dimension  cj_matrix_dim.   This  proc   } 

{  initialize  this  matrix  which  is  called  CJ_MATRIX.    } 

{  } 

procedure  INITIALIZE_CJ_MATRIX ; 

var 

j : integer; 

begin 

for  i:=  1  to  cj_matrix_dim  do 
begin 

for  j:=  1  to  cj_matrix_dim  do 
begin 

cj_matrix  [i,j]  :=  0; 
end; 
end ; 
end ; 
{******************************************************)- 

{  } 

{  INIT_REALIZATION_MATRIX                           } 

{  } 

{  In  this  program,  realizations  are  generated  one     } 
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{  interval  at  a  time  until  the  entire  realization  is   } 

{  constructed.   For  each  interval,  a  matrix  of  all    } 

{  known  realizations  is  read.   Additions  to  these     } 

{  partial  realizations  are  then  placed  in  a  different  } 

{  matrix.   This  proc  initialize  both  of  these        } 

{  matrices.  } 

{  } 

{**********+*******************************************} 

procedure  INIT_REALIZATION_MATRIX ; 
var 

i, 

j : integer; 

begin 

writeln  ('INITIALIZING  REALIZATION  MATRICES  '); 
writeln; 

for  1  :=  1  to  matrix_length  do 
begin 

for  j  :=  1  to  max_mtervals  do 
begin 

realization_x [l , j]  :=  0; 
realization_y [l , j]  :=  0; 
end; 
end; 

writeln  ('INITIALIZATION  COMPLETE'); 
writeln; 
end; 


{  } 

{  COPY  X_T0_Y  } 

{  } 

{  This  procedure  copies  the  contents  of  realization    } 

{  matrix  X  into  realization  matrix  Y.                 } 

{  } 

{*******************************************************}■ 

procedure  C0PY_X_T0_Y  (  x,y : integer) ; 
var 

i : integer; 

begin 

for   i    :=    1   to   max_mtervals   do 

realization.y [y , i]     :=   realization_x[x, l] ; 
end; 
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{  } 

{  COPY_Y_TO_X  } 

{  } 

{  This  procedure  copies  the  contents  of  realization    } 
{  matrix  Y  into  realization  matrix  X  } 

{  } 

{**********  +  *  +  *************  +  ****************************;- 

procedure  C0PY_Y_T0_X  (y,x: integer) ; 
var 

i : integer; 

begin 

for  i  :=  1  to  max .intervals  do 

realization_x[x, i]  :=  realization_y  [y ,i] ; 
end; 


{  } 

{   PRINT_MATRIX  } 

{  } 

{  Print_matrix  prints  out  the  most  recently  updated  } 

{  realization  matrix  (X  or  Y) .   If  the  flag  MATRIX  is  } 

{  set  to  1,  then  the  REALIZATION.X  was  the  last  } 

{  updated  and  is  subsequently  printed.   If  the  flag  } 

{  is  not  1,  the  Y  matrix  is  printed.   If  the  } 

{  diagnostic  flag  ARCHIVE  is  true,  output  goes  both  } 

{  to  the  screen  and  to  the  data  file  OUTFILEVAR.  } 

{  Else,  output  is  to  screen  only.  } 

{  } 

procedure  PRINT_MATRIX  (matrix,  length,  width: integer) ; 

var 

i, 

j : integer; 

begin 

for  i  :=  1  to  length  do 
begin 

writeln; 

if  archive  then 

writeln  (outf ilevar) ; 
for  j  :=  1  to  width  do 
begin 

if  matrix  =  1  then 
begin 

write  (realization_x [i , j] ) ; 
if  archive  then 
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write    (outfilevar, 

realization_x[i , j] ) ; 
end 
else 

begin 

write  (realization_y [i , j] ) ; 
if  archive  then 

write  (outfilevar, 

realizat ion_y [i , j] ) ; 
end; 
end; 
end; 
writeln; 
if  archive  then 

writeln  (outfilevar); 
end; 


{  } 

{  GENERATE.REALIZATIONS  } 

{  } 

{  This  procedure  produces  a  matrix  containing  all  } 

{  possible  realizations  given  an  initial  number  of  } 

{  interferers,  a  maximum  number  of  mterferers,  and  } 

{  the  total  number  of  intervals.   A  tree  method  is  } 

{  employed  to  build  up  a  realization  one  interval  at  } 

{  a  time.   The  first  interval  contains  the  initial  } 

{  number  of  interferers.   This  partial  realization  } 

{  is  placed  in  a  realization  matrix.   On  the  next  } 

{  pass,  the  partial  realization  is  read  and  the  next  } 

{  interval  is  added  to  it.   These  realizations  are  } 

{  then  placed  into  another  realization  matrix.   This  } 

{  process  continues  until  all  intervals  have  been  } 

{  filled  in  for  all  realizations.  } 

{  } 

procedure  GENERATE.REALIZATIONS  (mit.mterf  erers , 

max.interf erers , 

num_mtervals , 
min.mterf  erers  :  integer; 

var  active_matrix : integer ; 
var   active_matrix_length 
: integer) ; 

var 

temp, 
i, 

j> 
k, 
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x_list_length, 
y_list_length, 
interval_count : integer ; 

begin 

writeln  ('GENERATING  REALIZATIONS  FOR  ' ,max_interf erers , ' 
MAX  INTERFERERS  AND  ' ,num_ intervals , '  INTERVALS'); 
{      writeln  (init.mterf  erers,  '  '  ,max_interf  erers  ,  ' 
' ,num_intervals, '  ' ,min_interf erers) ;} 
writeln; 

if  archive  then 
begin 

writeln  (outf ilevar) ; 
writeln  (outfilevar, 

'REALIZATIONS  **************'); 
end; 

realization_x[l , 1]     :=    init.interf erers ; 
x_ list .length    :=    1; 
interval_count    :=    1; 
active_matrix    :=    1; 
active_matnx_length    :=    1; 

while   interval_count   <  num_intervals   do 
begin 

y_list_length    :=   0; 
{  writeln    ('within  x  to   y   loop');} 

for   i    :=    1   to   x_list_length  do 
begin 

if   realization_x[i , interval_count]    < 
max_mterf erers   then 
begin 

y_list_length    :=  y_list_length  +    1; 
copy_x_to_y    (l ,y_list_length) ; 

realization_y[y_list .length, interval_count+l] 

:=     realization_x[i , interval_count]+l ; 
end; 

if   realization_x[i ,interval_count]    > 
min.mterf erers   then 
begin 

y_list_length    :=  y_list_length  +    1; 
copy_x_to_y(i ,y_list_length) ; 

realization_y [y_list_length, interval_count+l] 

: =      realization_x [i , interval_count] -1 ; 
end; 
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end; 
interval_count  :=  interval.count  +  1; 
active_matrix  :=  2; 
active_matrix_length  :=  y_list_length; 

x_list_length  :=  0; 

if  interval_count  <  num_intervals  then 
begin 

for  k  :=  1  to  y_list_length  do 
begin 
{  writeln  ('  within  y  to  x  loop');} 

if  realization_y [k,interval_count] 

<  max_interf erers  then 
begin 

x_list_length  :=  x_list .length 

+  i; 
copy_y_to_x(k,x_list_length) ; 

real ization_x[x_list .length, interval _count+l] 

:=   realization_y [k, interval. count] +1 ; 
end; 
if  realization_y [k, interval_count] 
>  min_interf erers  then 
begin 

x_list_length  :=  x_list_length 

+  l; 

copy_y_to_x(k,x_list_length) ; 

realization_x[x_list_length , interval_count+l] 

: =   realization_y [k, interval.count] -1 ; 
end ; 
end; 
interval_count  :=  interval_count  +  1; 
active.matrix  :=  1; 

active_matrix_length  :=  x_list_length; 
end; 

end; 

writeln  ('There  Were  ' ,active_matrix_length, ' 
Realizations ' ) ; 
if  archive  then 
begin 

writeln  (outfilevar, 'There  Were 
' ,active_matrix_length, '  Realizations');  end; 

if  active_matnx  =  1  then 

print_matrix  (1 , active_matrix_length, 
interval_count) 
else 

print .matrix  ( 2, act ive_matrix_ length, 
interval_count) ; 
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end; 


{  } 

{  FILTER  } 

{  } 

{  This  procedure  accepts  a  realization  matrix  which  } 

{  contains  all  possible  realizations  given  a  specific  } 

{  combination  of  interferers.   Each  realization  in  } 

{  the  matrix  is  checked  for  the  invalid  condition  of  } 

{  overcount.   An  overcount  occurs  when  a  realization  } 

{  has  indicated  that  a  certain  number  of  interferers  } 

{  has  departed  but  yet  somewhere  in  the  realization  } 

{  the  interferer  count  goes  beyond  the  count  of  the  } 

{  remaining  interferers.   Realizations  which  are  } 

{  filtered  out  are  discarded.   All  good  realizations  } 

{  are  placed  into  the  currently  inactive  realization  } 

{  matrix.  } 

{  } 

procedure  FILTER  (var  length : integer ;  intervals, 
interferers : integer ; 

var  active_matrix: integer) ; 

var 
i, 
J. 

P» 
max , 

num_f mal_interf erers : integer; 

invalid_realization : boolean; 

begin 

writeln; 

writeln  ('FILTERING  THE  REALIZATION  MATRIX  '); 

writeln; 


J  :=  1; 

if  active_matrix  =  1  then 
begin 

active_matrix  :=  2; 
for  i  :=  1  to  length  do 
begin 

invalid_realization  :=  false; 
num_f mal_interf erers  := 

realization_x[i , intervals] ; 
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if  num_f inal_mterf erers   =   interferers  then 
begin 

if   intervals   >   2  then 
begin 

max    :=  max_interf erers; 
for  p    :=   1   to    (intervals   -    1)    do 
begin 

if   realization_x[i ,p]    > 

realization_x[i,p+l]    then 
max    :  =  max   -    1 ; 
if   realization_x[i,p+l]    >  max 
then 

invalid_realization 
:=  true; 
end; 
end; 
if    invalid_realization   =   false  then 
begin 

copy_x_to_y    (i,j); 
{  print _matrix(2 , j ) ; } 

j    :=  j    +   1; 
end; 
end; 
end; 
end 
else 

begin 

active_matrix    :=    1; 
for   l    :=    1    to    length   do 
begin 

mvalid_realization    :=  false; 
num_final_mterf erers    :  = 

realization_y[i , intervals] ; 
if   num_fmal_interf erers   =    interferers   then 
begin 

if   intervals   >   2  then 
begin 

max    :=   max_mterf erers ; 
for  p    :=    1   to    (intervals   -    1)    do 
begin 

if   realization_y [l ,p]    > 
realization.y [i,p+l]    then 
max    :  =  max   -    1 ; 
if   realization_y [i,p+l]    >  max 
then 

invalid_realization 
:=  true; 
end; 
end; 
if    mvalid.realization  =   false  then 
begin 


copy_y_to_x  (i,j); 
{  print_matrix  (l,j);} 

j  :-  j  +  1; 

end; 
end; 
end; 
end; 
writeln; 

writeln  ('THERE  WERE  ',j-l,'  VALID  REALIZATIONS  FOUND'); 
writeln; 
if  archive  then 
begin 

writeln  (outf ilevar) ; 

writeln  (outf ilevar , 'FILTER  ********************'); 
writeln  (outf ilevar , 'There  Were  ',j-l,'  Valid 
Realizations ' ) ; 
end; 
length  :=  j-1; 

print.matrix  (act ive_matrix, length, intervals) ; 
writeln; 
if  archive  then 

writeln  (outf ilevar) ; 
end; 


{  } 

{  GET.Cj  } 

{  } 

{  GET_Cj  reads  in  all  valid  realizations  individually  } 

{  and  determines  the  maximum  number  of  interferers  } 

{  in  each  realization.   A  tally  is  kept  of  how  many  } 

{  realizations  have  a  given  maximum  number  of  } 

{  interferers.  } 

{  } 

procedure  get_cj  (active_matrix,  active_matrix_length, 

n :  integer) ; 

var 
i, 

j. 

max : integer; 

begin 

writeln  ('FINDING  Cjs  FOR  A  REALIZATION  WITH  ' ,n, ' 
PACKETS ' ) ; 

if  archive  then 
begin 
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writeln  (outf ilevar) ; 

writeln  (outf ilevar , 'CJ  CALCULATION  *********'); 

writeln  (outf ilevar) ; 
end; 

for  i  :=  1  to  active_matrix_length  do 
begin 

if  active_matrix  =  1  then 
begin 

max  :=  realization_x[i, 1] ; 
for  j  :=  2  to  max_intervals  do 
begin 

if  max  <  realization_x[i , j]  then 
max  :=  realization_x[i , j] ; 
end; 
cj.matnx  [n,max]  :=  cj_matrix  [n,max]  +  1; 
if  archive  then 
begin 

writeln  (outf ilevar) ; 

writeln  (outfilevar,  'MAX  FOR  THIS  PASS 

IS  '.max); 
writeln  (outf ilevar , 'CUMULATIVE  MAX  FOR 
I=',max,'  IS  ' ,cj .matrix [n, max] ) ; 
end; 
end 
else 

begin 

max  :=  realizat ion_y  [l , l] ; 
for  j  :=  2  to  max.mtervals  do 
begin 

if  max  <  realization.y [i , j]  then 
max  :=  realization_y [l , j] ; 
end; 
cj.matrix [n,max]  :=  cj .matrix [n, max]  +  1; 
if  archive  then 
begin 

writeln  (outfilevar) ; 

writeln  (outfilevar,  'MAX  FOR  THIS  PASS 

IS  ' ,max) ; 
writeln  (outf ilevar , 'CUMULATIVE  MAX  FOR 
I=',max,'  PKTS=',n,'  IS 
' , cj .matrix [n, max] ) ; 
end; 
end; 
end; 
end; 


{  } 

{   PRINT.CJ  } 
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{  } 

{  This  procedure  prints  out  the  matrix  of  Cj  values   } 
{  constructed  by  the  procedure  GET_Cj .  } 

{  } 

{******************************************************} 

procedure  PRINT.CJ; 

var 

i, 

j : integer; 

begin 

writeln  ('     I   123456789   10'); 
for  i:=  1  to  cj_matrix_dim  do 
begin 

write  (i,'    I ') ; 
if  archive  then 

write  (outf ilevar , l , '    I'); 
for  j :=  1  to  cj_matrix_dim  do 
begin 

write('  ' ,cj_matrix[i , j] , '  '); 
if  archive  then 

write  (outf ilevar , '  '  ,cj_matrix[i ,  j]  )  ; 
end ; 
writeln; 
if  archive  then 

writeln  (outf ilevar) ; 
end ; 
end; 


-[  ++************+***+**  *  *  *  *  *  ***********: *  *************  *  *  +;■ 
{  } 

{  COMB  } 

{  } 

{  Given  a  total  number  of  mterferers,  every  } 

{  combination  of  mterferer  packet  type  must  be  } 

{  determined.   This  procedure  accepts  the  total  } 

{  number  of  mterferers  n  and  generates  every  } 

{  combination  of  the  four  interferer  types  which  } 

{  add  up  to  the  total  n.   This  is  basically  an  } 

{  enumeration  of  the  number  of  combinations  of  4  } 

{  things  taken  n  at  a  time.   This  algorithm  is  } 

{  derived  from  Reference  7.  } 

{  } 

-T***********  **************  *****************************}- 

procedure  COMB  (n,k : integer ;  var  r:index;  var  mtc :boolean) ; 
label   10,  11,  15,  20,  30,  31; 
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var 
i, 

J. 

w, 
x, 

y> 

z, 

1 : integer; 

begin 


10: 

if 

mt 

c   =   true   then 

begin 

goto   20; 

en 

d 

el 

se 

be 

igin 

r[l]     :=   n; 
t    :=   n; 
h    :=   0; 

11: 

15: 

20: 
30: 

for   1    :=   2   to   k   do 

r[l]     :=   0; 
if   r[k]    =   n   then 

mtc    :=   false 
else 

mtc    :=   true; 
goto   31 ; 
if   t    >    1   then 

begin 

h    :=   0; 

end; 
h    :=   h   +    1; 
t    :=   r[h]; 
r[h]     :=   0; 
r[l]     :=   t    -    1; 
r[h+l]     :=   r[h+l]    +    1; 
goto    15; 

31: 

en 
1 

d; 
:=   1; 

end; 


{  } 

{  GET.COMBINATIONS  } 

{  } 

{  This  procedure  controls  the  execution  of  the  } 

{  procedure  COMB  and  generates  the  values  w,  x,  y,  } 

{  and  z.   These  variables  correspond  to  early-late,  } 

{  early-early,  late-early,  and  late-late  packets  } 

{  respectively.   All  combinations  given  a  total  } 
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{  number  of  interferers  n  are  placed  into  the  matrix  } 
{  vector.  } 

{  } 

procedure  GET_COMBINATIONS  (  var  vector :comb_matrix; 

n: integer;  var  vector.length 
integer) ; 

var 
i, 

J. 
w, 
x, 

y> 

z: integer; 

mtc :boolean; 

v: index; 

begin 

for  1  :=  1  to  4  do 
v[i]  :=  0; 

i  :=  0; 

mtc  :=  true; 

if  archive  then 

writeln  (outf llevar, 'COMBINATIONS 

*** *  *  * *  * *  *  ******  *  *  +  +  >  )  ■ 

while  mtc  do 
begin 

if  l  =  0  then 

mtc    :=  false; 
i    :=    i   +    1; 
comb    (n,4, v,mtc) ; 
for  j    :=    1    to   4   do 

vector[i,j]     :=   v[j]  ; 


w 

=   v[l] 

X 

=   v[2] 

y 

=   v[3] 

z 

=   v[4] 

writeln  (w,J  ' ,x, '    ' ,y , '    \z);} 
if  archive  then 
begin 

writeln  (outf ilevar ,w, '  ',x,'  ' ,y , '    ' ,z) ; 
end; 
end; 
vector_length  :=  i; 
if  archive  then 
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writeln 
(outfilevar, '*********************************') ; 
end; 

{  } 

{  MAIN  } 

{  } 

{  Program  execution  begins  with  a  query  for  an  } 

{  archive  data  file.   If  the  archive  option  is  chosen  } 

{  all  screen  I/O  will  be  echoed  to  the  file  given.  } 

{  The  next  program  query  is  for  the  maximum  number  of  } 

{  packets,  n.   The  program  will  then  generate  all  } 

{  combinations  of  interferer  packet  types  for  a  given  } 

{  value  of  n.   Each  combination  is  then  fed  one  by  } 

{  one  into  the  GENERATE.REALIZATION  procedure  where  } 

{  all  possible  realizations  are  generated.   This  list  } 

{  of  realizations  is  then  filtered  by  FILTER  to  cull  } 

{  out  invalid  realizations.   A  matrix  of  valid  } 

{  realizations  is  fed  to  GET_Cj  in  order  to  generate  } 

{  the  Cj  table.   This  process  is  repeated  for  all  } 

{  values  of  n  between  1  and  the  maximum  given  by  the  } 

{  user.  } 

{  } 


begin 

archive  :=  false; 

write  ('DO  YOU  WANT  A  DATA  ARCHIVE  (y  or  n)  ->  >) ; 
readln  (ans) ; 
if  ans  =  'y'  then 
begin 

archive  :=  true; 

write  ('ENTER  DATA  FILE  NAME  ->  '); 

readln  (f ile_name) ; 

assign  (outf ilevar,  f ile_name) ; 

rewrite  (outf ilevar) ; 
end; 

write  ('ENTER  THE  MAX  #  OF  PKTS  PER  REALIZATION  ->  '); 
readln  (runs) ; 

initialize_cj .matrix ; 
for  n  :=  1  to  runs  do 
begin 
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for  i  :=  1  to  220  do 
for  j  :=  1  to  4  do 
vector[i,j]  :=  0; 

if  archive  then 
begin 

writeln 
(outf ilevar , '*********************** 

*  +  *  *  *  +  +  +  *******  +  **********+'  )  - 

writeln    (outf ilevar , 'NUMBER  OF  PACKETS  =    ' ,n) ; 
writeln; 
end; 

get_combmations    (vector  ,n,vector_length)  ; 
for   i    :=    1   to  vector_length  do 
begin 

init_interferers    :=  vector [i, l]+vector[i ,2] ; 
f inal_interf erers    :=  vector [i , l]+vector[i ,4] ; 
max_interf erers    :=   vector [i,l]    +  vector [i, 2]    + 

vector[i,3]    +  vector[i,4]; 
min_interf erers    :=   vector[i,l]; 
num_intervals    :=  vector[i,2]    +   2*vector [i ,3]    + 
vector  [i ,4]    +    1; 


writeln 

(outf 

writeln 

(outf 

writeln 

(outf 

writeln 

(outf 

writeln 

(outf 

writeln 

> 

(outf 

writeln 

(outf 

writeln 

(outf 

writeln 

(outf 

writeln 

(outf 

writeln 

(outf 

writeln 

(outf 

writeln 

(outf 

if  archive  then 
begin 
ilevar) ; 
ilevar) ; 

ilevar,  ' ' ) ; 

ilevar) ; 

ilevar, 'PACKET  CONFIGURATION  IS  ', vector [l , 1] , 

'  ' ,vector[i,2] , ;  ', vector [i, 3] , ; 
,  vector  [l  ,4]  )  ; 
ilevar) ; 

ilevar, 'INITIAL  PKTS  = 
ilevar, 'FINAL  PKTS 
ilevar, 'MAX  #  PKTS 
ilevar, 'MIN  #  PKTS 
ilevar, 'NUM  INTERVALS= 
ilevar) ; 

ilevar, ' 

end; 


' , init_interf erers)  ; 
' ,final_interf erers) ; 
' ,max_interferers) ; 
' ,min_interf erers) ; 
' ,num_intervals) ; 

>). 


mit_realization_matrix; 

generate_realizations    (mit_interf erers , 
max_interf erers , 

num.intervals,    min_interf erers , 
active_matrix, 
active_matrix_length) ; 
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filter  (active_matrix_length,  num_ intervals , 
f inal.interf erers , 

active_matrix) ; 

get_cj  (active_matrix,  active_matrix_length, 
max_interf erers) ; 
end; 
end; 
print_cj ; 
close  (outf ilevar) ; 
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APPENDIX  B 

DERIVATION  OF  THROUGHPUT  EQUATION 
FOR  VARIABLE  PACKET  ALOHA 

This  appendix  presents  the  derivation  of  the  throughput  equation  for  variable 
packet  Aloha  with  power  capture. 

As  seen  in  Chapter  II  equation  (2.9).  the  probability  of  a  tagged  packet  success- 
fully being  transmitted  is 

n     N 


p.(*,«)  =  e-^'ll  +  f;/(n)2-<"+"lg(x  +  ")]"  £<?>)]  (B.l) 

Assign  the  following  identity 


fc(n)  =  /(n)2-t»+1^X;Ci(n)  (B.2) 


ill 

J  =  0 


As  described  in  Chapter  II  equations  (2.10)  and  (2.11).  the  conditional  channel 
departure  rate  and  the  channel  throughput  are  given  by  the  following: 

d>(x,u)  =  g(x  +  u)pa(x,u)  (B.3) 

/■oo      roc 

S  =  —  I      /      (f>(x,u)ci(x)a(u)dxdu  (B.4) 

2  Jo     Jo 

Utilizing  the  binomial  identity 


(*  +  «)"  =  £      I   )xkun-k  (B.5) 


4C 


and  expanding  the  exponential  term  in  (B.l),  the  equation  for  ps(x,u)  can  be  rewrit- 
ten as 


ps{2\u)  =  e-9Xe-gu 
Expansion  of  4>(x,u)  in  terms  of  ps(x,u)  results  in 


1  +  £  Kn)  E 


n      xku-k 


(B.6) 


<f>(x,  u)  =  g(x  +  u)e  gxe  9U 
Multiplying  produces 


i  +  E^WE 

n-\  k=0 


(B.7) 


0(x,u)     =     [^ze^e-5*  +  gue-9Ue~9X] 
+[gxe-9Ue~9X  +  gue-9Ue-ax] 


EMn)E 


n  =  l 


A=0 


X   U 


(B.8) 


Using  this  definition  of  (f>(x,  u)  in  the  equation  for  S  (B.4)  produces  two  separate 
double  integrals.  The  first  is  equivalent  to  the  equation  (27)  in  Ref.  2.  Assuming 
that  Pi(x)  =  0.  pziv)  —  0  and  that  the  density  function  a(x)  is  exponential,  this 
first  term  for  S  is  found  in  Ref.  2  to  be 


5i  = 


G 


(B-9) 


vG  +  1)3 

The  remaining  double  integral  remains  to  be  evaluated.  This  term  represents  the 
increase  in  throughput  due  to  power  capture.  The  term  is  written  here  as 


roc      roc 

S2     =     -  /       /     [gxe-9Te-9U  +  gae^e'9*] 
2  Jo     Jo 


±h(n)±(l)x^ 

n=i  k=0  V    k    J 


a{x)a(u)dxdu 


(B.10) 


Manipulating  terms  in  the  above  equation  produces 


52    =    - 


[      roc-      roc 


ax    —  o 

xe  y  e  - 


"E 

fc=0  L 


72     1   „*„«-* 


TOO       TCxD  n 


Jt=0  L 


n     \      kn-k 
k         XU 


a(x)a(u)dxdu 
a(x)a(u)dxdu>  (B.ll) 


Further  collection  of  like  terms  produces 


1     °°  (  n    /  n  \    r°°  r°° 

S2    =     -gY.h(n)\12  [I       /     xk+1e-9Xa{x)dx-         un-ke-gua{u)du 

+  J2  (  nk   )  r un-k+1e-9Ua(u)du-  f°°  xke-9Ta{x)dx\  (B.12) 

Each  of  the  individual  integrals  above  can  be  solved  using  the  following  identity 
[Ref.  5:p.  310]. 


roc 

L  xn 

Jo 


e'^'dx  =  nhr{n+l) 


B.13) 


Additionally,  the  following  substitutions  can  be  made  because  packet  size  is  as- 
sumed to  be  exponentially  distributed  and  because  the  quantity  G  is  given  in  pack- 
ets/average packet  length. 


a(x)  =  —  e 
x 


c/r 


;b.h) 


G  =  gx 


(B.15) 


Combining  the  integral  identity  in  equation  (B.13)  with  the  definition  of  a(x) 
produces  this  solution  for  the  first  integral  term  in  equation  (B.12). 

/•OO  1         roo  . 

/     xk+1e-9Xa{x)dx    =     -         xk+le^9+  =  )xdx 
Jo  x  Jo 
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2:  X 

(fc  +  1)!  gx  +  1     (k+2) 

X  X 

=    !^(C+l)-<t+2>  (B.16) 

The  second  term  in  equation  (B.1'2)  is  evaluated  in  the  same  manner  as  above. 

/      un-ke-gua{u)du     =     -  /      un-ke-(9+  =  )udu 
Jo  u  Jo 

_     (w  -  fc)!       |  j.     (n_fc+1) 

=     (lzJQi(^±i)-(n-*+i)  (B.17) 

u  u 

In  this  model,  all  packets  come  from  the  same  distribution,  therefore  x  =  77.   This 
modifies  the  above  equation. 

r  un-ke-9Ua(u)du  =  (l'_r_AM!(G  +  ir(7'^+1)  (B.18) 

The  third  integral  term  in  equation  (B.1'2)  is  reduced  as  follows 

jTV-*«e-»«a(u)du  =  i|^±±ii!(G+  l)-<"-1+2'  (B.19) 

The  final  integral  term  is  as  follows 

I"  xne~axa{x)dx  =  4t(6'  +  l)~{k+l)  (B.20) 

Jo  X 

The  total  throughput  of  the  system  is 

S  =  S\  +  S2  (B.21) 

Collecting  equations  (B.9),  (B.ll).  and  (B.16)-(B.20)  produces 
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s   = 


G 


+  (^,1,)'(g+ i)-"-w»|r(g+ i)-'t+1' 


^iT(G  +  l)-,fc+2) 


This  simplifies  to 


5    = 


(GTTpn9SAWter+1)i£^r- 
+-'(-*+D^±s^]} 


Further  simplification  produces 


(B.22) 


(B.23) 


5 


G 


+  ^  E  *0 


(G+1)_(n+3>     n 


\3         9 


-(n+1) 


£n!(n  +  2) 


(B.24) 


(G  +  l,  n=1  fc=0 

The  final  summation  in  the  above  equation  is  over  the  range  of  k  —  0  to  oc. 

Because  there  are  no  k  terms  within  the  summation,  the  entire  summation  reduces 
to  n\n(n  +  2).  Combining  terms  produces 


S 


G         ,    1      ~   /(n)   "  ^(C;+1)-(n+3) 


n  =  l 


j  =  0 


Combining  terms  and  utilizing  the  identity  G  =  gx  produces  the  following  solu- 
tion for  the  throughput  of  variable  packet  Aloha  with  power  capture. 


V0        '        -1     '  71  =  1  ~  j=0 


(B.26) 
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APPENDIX  C 

DERIVATION  OF  THROUGHPUT  EQUATION 
FOR  SELECTIVE-REPEAT  ALOHA 

This  appendix  presents  the  derivation  of  the  throughput  equation  for  the  variable 
packet  Aloha  network  with  power  capture  and  selective-repeat  data  link  control 
presented  in  Chapter  III. 

From  Chapter  III  equation  (3.2),  the  probability  of  tagged  packet  success  to  be 


P,(Z)    =    £-»(»"  +  '"> 


l  +  f/(n)2-("">^ +  "*)]"£  Ci(, 


n  =  l 


n'. 


J=0 


(C.27 


The  conditional  probability  of  successful  minipacket  transmission  is  given  in  equa- 
tion (3.3)  of  Chapter  III  to  be 


q{n,y)  =  ps{z)[l  -  Pi{m)][l  -  p2{y)\ 


(C.28) 


Also  given  in  Chapter  III  is  the  conditional  minipacket  departure  rate  in  the  interval 
of  n  +  1  minipackets.  This  equation  is  rewritten  here. 


<f>m{z,y)  =  mg{z  +  l)q{z,y) 


(C.29) 


Equation  (3.5)  in  Chapter  III  gives  the  expression  for  the  average  minipacket  through- 
put and  is  rewritten  here. 


/•oo 

Sm  =  =— r  /     <f>m(z,y)dy 
z  +  1  Jo 


(C.30) 
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The  equation  for  minipacket  throughput  given  above  is  now  expanded  to  produce 
the  following: 


^  —  mzg 


z  r°°  .  , 

Sm    =    — -mg[l  -  Pl(m)]e-ms        [1  -  p2(y)]b(y)dy  ■  J>  +  l)D(z) 
2  +  1  Jo  f^[ 

—  -oo  CO         oo 

+  —-rmg[l  -  Pl(m)]e-^  /     [1  -  p2(y)}b(y)dy  ■  £  •  £  /(n)2"<n+1> 

>(m*"  +  m)1"  £  [C»]  (,  +  i)2>Me-*  (C.31) 

where  D{z)  is  the  probability  that  a  packet  consists  of  z  minipackets. 

The  first  term  in  the  above  equation  is  the  same  as  the  equation  evaluated  in 
Ref.  2  equation  (36).  The  resulting  throughput  equation  for  this  term  is  given  in 
Ref.  2  equation  (37b).  Call  this  part  of  the  total  throughput  S\.  Then  from  Ref.  2, 


_  ffC,y(2-E)[l-p,(m)] 
-''  "      5(5  +  l)u-'(l-«)»  (C'32) 

With  the  following  definitions: 


(C.33) 
(C.34) 

(C.35) 
(C.36) 

Let  52  be  the  remaining  term  of  the  total  throughput  given  in  equation  (C.31). 
By  separating  terms  and  noticing  that  [g(mz  +  m)]n  =  (gm)n(z  +  l)n,  S2  is  rewritten 
as 


c  mzg  n  /     m  -ma 

S2     =     — — [1  -pi(m)je      V 

z  -\- 1 


52 


1 

u>    =     1  -  - 

z 

G 

v    =    we   * 

G       =       7722^ 

/•CO 

a     =      Ml- 

./0 

-  P2{y)]b{y)dy 

n=l 


■mzg 


(C.37) 


(nm)n    N  °° 

/(n)2-(n+i)l^L  £  [Cj{n)]  .  J2(z  +  lT^D(z)e 

As  stated  in  Chapter  III  equation  (3.7),  assume  that  D(z)  is  geometrically  dis- 
tributed and  is  defined  as  follows: 


D(z)=  =(l-z)2_1  (C.38) 

z  z 

Utilizing  the  equations  for  G,  D(z),  and  w  the  equation  for  52  can  be  written  as 


follows: 


s, 


G 


3+1 
■i)/(n)2 


[l-Pl(777)]e-(fV 


/  G_  \  n     A' 
9-("+l)ixi 


,|G- 


n=l 


71 ! 


Eto^EC^1)^1^"1'"^     (C39) 


j  =  0 


2  =  1 


Assign  the  following  identity: 


Hn)  =  f(n)2-{n+1)^^IECj(n)  (C40) 

At  this  point,  consolidate  the  exponential  terms  and  divide  out  two  of  the  w 
terms  to  equalize  the  exponents  within  the  summation  over  z. 


G 


S2  =  ^[l-Pl(m)]aj:h(n)w-2Y,(z  +  l)n+1™Z+1e-mz+1] 

Z  +    ]  n  =  l  z=0 


Assign  the  following  identities: 


(C.41) 


=    (*  +  l) 
=     (n  +  1) 


(C.42) 
(C.43) 


This  produces  the  following 
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&  = 


G 


+  1 


[l-ft(ro)]a£>(n)uT2][>V 


(C.44) 


n=l  k=2 

At  this  point,  modify  the  summation  over  k  such  that  the  range  is  from  A;  =  1 
to  oo.  To  do  this,  subtract  off  from  the  total  of  the  summation  the  value  of  the 
summation  for  k  =  1  which  is  v  in  this  case.  Therefore,  the  above  equation  is 
changed  to 


s2  =  - 


G 


r  +  1 


[1  -pi(m)]a^2  h(n)w~ 


n  =  l 


Ufc=i 


t..k 


(C.45) 


Given  that  \v\  <  1,  use  the  following  mathematical  identity  to  solve  the  above 
equation  [Ref.  6:  p.  142]. 


r- 1 


(r-ty 


v-(1--^('+1,^1^^^^-,+l)! 


(C.46) 


k=l 


Using  the  above  identity,  the  final  form  of  52  is  given  as 


•S'2    = 


G 


z  +  \ 


[1  -  Pi{m)}aJ2  h(n)w 


71  =  1 


V  -  k)> 


(l-)-^(«+l)l^l(-l)*H(«_4+I)1 


—  V 


(C47) 


The  total  throughput  for  selective-repeat  Aloha  with  power  capture  is  given  as 


S  -  S1  +  52 


(CAS) 
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APPENDIX  D 

ANALYTICAL  AND  NUMERICAL  RESULTS  FOR 
NEAR/FAR  EFFECT  CAPTURE 

This  appendix  provides  the  details  of  the  analytical  and  numerical  methods  used 
to  determine  the  values  for  Pr{ capture |/z}  in  Chapter  III.  First,  the  derivation  of 
Pr{capture\u  =  1}  is  provided.  This  is  followed  by  a  MATHCAD  file  which  was 
used  to  calculate  capture  for  n  =  2  and  3. 

A.  ANALYTICAL  METHOD 

The  case  of  v  =  1  is  sufficiently  simple  to  conduct  a  straightforward  mathematical 
derivation  of  the  capture  probability  given  the  near/far  effect.  From  Chapter  IV  the 
following  is  known: 


Z=y  (D.49) 


/A'(-r)  =  ^^(2/a+1)  (D.50) 


fv(y\n)  = 


aGy 


,  I, 


(D.51) 


For  the  case  of  n  =  1,  both  y  and  x  are  defined  between  the  values  of  r"^  and 
oo.  Also  from  Chapter  IV,  equation  (4.10),  the  definition  of  fz(z\n).  From  this 
equation,  it  is  clear  that  both  yz  and  y  must  be  greater  than  r~£r  given  the  valid 
ranges  for  fxix)  an<J  /v(y|7')-  This  implies  that  for  values  of  z  greater  than  1,  the 
valid  range  for  y  is  y  >  r~"x.  Therefore,  the  equation  for  fz{z\n)  with  the  proper 
limits  of  integration  is 
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fz(z\n  =  1)     =      /       yfx{yz)fv(y\n  =  l)dy 


-l 


max 
OO 


»^(y*)-w*+,,^»-("*+,)* 


VaG/  A-s, 


a 


G 


Remembering  that  G  =  Kr2max  gives  us 


fz(z\n  =  1 


ar4 

max 


.-(2/or  +  l)   „-4/ 


.'/ 


ir-(2/a  +  l) 
Q 


(D.52) 
(D.53) 
(D.54) 
(D.55) 


(D.56) 
(D.57) 


To  determine  the  capture  probability  from  the  above  distribution,  simply  inte- 
grate from  the  capture  threshold  70  to  infinity. 


Pr{capture\n  —  1}     = 


Pr{z>l0\n  =  1} 

(D.58) 

/     -z-(2l°^dz 

(D.59) 

1      -2/q 

^7o 

(D.60) 
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B. 


NUMERICAL  METHODS 


This  document  is  a  MATHCAD  file  used  to  numerically  calculate 
the  values  of  P{capture|n}  for  values  of  n  from  2  to  3 .   The 
general  approach  is  to  use  a  numerical  convolution  to  determine 
the  n'th  convolution  of  f (y) .   These  convolution  results  are 
then  used  in  numerical  integrations  to  generate  the  P{capture|n} 
results. 


Initialize  and  Define  Variables 


M  =  1 

<x  S  4 

(  M  corresponds  to  r   ) 

max 

2 

G  a  t  M 

D  =  4 

(  2  <  D  <  32  sets  the  nu 

i  =  0  .  . 

32 

D  -  1 

2  i 

+ 

1 

V   — 

i     2 

D 

f(t) 


2 
-+1 


=  2 t       •  J*  (10  -  t)  -  *  [M    -  tj. 


<x  G 
h(t)  =    f(t) 
P(t)  =  f(t) 


f(t)  corresponds  to  f(w)  in  the  model 


The  following  eguation  performs  a  discrete  convolution  on  the 
functions  defined  as  X  and  Y. 


conv(X,Y)  = 


icfft 


[\  16-  D- 


(cfft(X)  cfft(Y)) 


a 
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Step  1.    Convolve  f(t)  with  itself  to  get  f(y|n=2)  (call  this 
C)  and  plot  the  resulting  functions. 


C  =  conv 


|h(v),p(v)J 


0.06 


-0.01 


40 


Step  2.    Convolve  f(y|n=2)  with  f(t)  to  get  f(y|n=3)  (called 
D)  and  plot  the  relulting  function. 


D  =  conv 


l_h(v),c| 


0.02 


D 


-0.01 


V 


40 


Step  3 .    Continue  in  the  above  manner  for  as  many  cases  of 
f(yln)  as  needed. 
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Step  4.    Calculate  f(z|n)  using  a  linear  interpolation  of 
the  numerical  value  for  f(y|n). 


2 

- 

-+1 

IT 

ex 

f(z)     :=   2 

Z 

ex  G 


.  30 


<y. 


• linterp(v,C,y)  dy 


-ex 


2M 


Step  5.    Calculate  the  capture  probability  given  n  known 
interferers  by  integrating  f(z|n)  (f(z)  in  this 
example)  from  the  capture  threshold  (3  in  this 
example)  to  a  maximum  useful 1  limit  which  is 
determined  by  the  range  of  values  for  which  f(y|n) 
is  valid.   In  this  case,  30  is  the  upper  limit. 


■-1 


30 


f(z)  dz 


Z  =  0.057 


Step  6.    Continue  on  in  the  same  manner  as  above  to 

calculate  the  capture  probability  for  different 
values  of  n.   The  following  equation  defines  f(z|n) 
for  the  case  where  n=3. 


9(h) 


ir 

2 h 

ex   G 


2 

-+1 
ex 


r»30 


ex: 


• linterp(v,D,y)    dy 


-ex 


3-M 


30 
3 
H  =   0.022 


•  30 

H  :=      g(h)  dh   This  equation  calculates  P{ capture | n=3 } 
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